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Solid CS 2 is superficially similar to CO 2 , with the same Cmca molecular crystal structure at low 
pressures, which has suggested similar phases also at high pressures. We carried out an extensive first 
principles evolutionary search in order to identify the zero temperature lowest enthalpy structures of 
CS 2 for increasing pressure up to 200 GPa. Surprisingly, the molecular Cmca phase does not evolve 
into /3-cristobalite as in CO 2 , but transforms instead into phases HP2 and HP1, both recently 
described in high pressure SiS 2 - HP1 in particular, with a wide stability range, is a layered P2i/c 
structure characterized by pairs of edge-sharing tetrahedra, and theoretically more robust than all 
other CS 2 phases discussed so far. Its predicted Raman spectrum and pair correlation function agree 
with experiment better than those of /3-cristobalite, and further differences are predicted between 
their respective IR spectra. The band gap of HPI-CS 2 is calculated to close under pressure yielding 
an insulator-metal transition near 50 GPa in agreement with experimental observations. However, 
the metallic density of states remains modest above this pressure, suggesting a different origin for 
the reported superconductivity. 


I. INTRODUCTION 


The crystal structure of even extremely stable molec¬ 
ular carbon compounds like CO 2 is known to transform 
radically at high pressures where, above ~50GPa, the 
carbon coordination is found, both experimentally 1 and 
theoretically 2 to switch from two in the molecular struc¬ 
ture Cmca to four in the /3-cristobalite structure.-® Sim¬ 
ilar pressure-induced structural transformations can rea¬ 
sonably be expected to occur in a compound such as 
CS 2 , which is the focus of the present study, and which 
presents obvious similarities to C0 2 Indeed CS 2 adopts 
at zero pressure the same Cmca molecular crystal struc¬ 
ture as CO 2 at moderate pressures between 15 and 50 
GPa. However, the binding in CS 2 is much weaker 
than in CO 2 , with a smaller electronic band gap and a 
positive formation enthalpy of about 88.7kJ/mol as op¬ 
posed to a large negative one of —393.509kJ/mole for 
CO 2. 56 Thermodynamically, that makes decomposition 
and phase separation into elementary carbon and sul¬ 
fur a thermodynamic necessity for CS 2 at equilibrium al¬ 
ready at ambient pressure and presumably even more so 
at higher pressures. In spite of that intrinsic thermody¬ 
namical metastability, solid CS 2 phases do exist, clearly 
for kinetic reasons, and are reported at ambient pres¬ 
sure not to dec ompose in measurable times, at least be¬ 
low « 560 Kj® a temperature rising even further at high 
pressures, possibly up to 1000 K at 70 GPaf' At high pres¬ 
sures but low temperatures, X-ray data have shown evi¬ 
dence of structural transitions of CS 2 from the molecular 
phase to polymeric phases with C-S coordinations rising 
from two (Cmca) to three (CS 3 ) to four (CS^Pand fur¬ 


ther. With the exception of the CS 3 phase near 20 GPa, 
the high pressure structural behavior has been postulated 
so far to be similar to that of CO 2 , implying (not unrea¬ 
sonably) that CS 4 could be /3-cristobalite. In CS 2 a de¬ 
tailed interpretation of high pressure experimental data 
is further complicated by a large amount of structural 
disorder, particularly in the CS 4 phase where only broad 
rather than sharp Bragg peaks are present in the diffrac¬ 
tion pattern. ^ With that rationalization of the fourfold 
coordinated state of CS 2 near 50 GPa through simple 
analogy with CO 2 , density functional theory (DFT) stud¬ 
ies have considered /3-cristobalite (I42d, also referred to 
as Chalcopyrite) and tridymite (P2i2i2i), for compar¬ 
ison with experimental data. Lacking sharp diffraction 
peaks or other distinctive features, that comparison nev¬ 
ertheless appears somewhat elusive. The tridymite struc¬ 
ture shows slightly better agreement with experiment 
but represents a thermodynamically less likely candidate 
than /3-cristobalite since its calculated DFT enthalpy is 
0.3-0.4eV/molecule higher.' The electronic structure of 
either fourfold coordinated crystal structure agrees with 
the observed metalization in the region from 40 GPa up¬ 
wards !- 7 9 Besides that, in more recent experiments per¬ 
formed on the high pressure metallic phases of CS 2 , re¬ 
sistivity showed evidence of superconductivity at 4-6 K 
over a broad pressure range from 50 to 172GPaP 

From the theoretical point of view it is of course inade¬ 
quate to merely trust the analogy between CS 2 and CO 2 
and extend it even to very high pressures. To expand 
somewhat our view we instead consider that CO 2 and 
CS 2 are members of a broader family of IV-VI AB 2 com¬ 
pounds including SiC> 2 , silica, and SiS 2 , neither of which 
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is molecular at zero pressure. Silica is very well known 
for a number of tetrahedrally coordinated polymorphs 
including a-quartz, a- and /3-cristobalite, tridymite, co- 
esite, etc. Much less studied until recently, SiS 2 ex¬ 
hibits totally different phases starting off at zero pres¬ 
sure with an orthorhombic (so-called NP) crystal struc¬ 
ture made up of edge-sharing carbon-centered tetrahedra 
(see Ref flUj) Very recently the high pressure phases of 
SiS 2 were experimentally characterized 10 and a whole se¬ 
quence of phases HP1, HP2, HP3, were described where 
the tetrahedra did not disappear but simply changed 
their mutual connectivity. 

Based on the limited available high-pressure data, solid 
CS 2 is at the moment the least understood member of 
this highly important class of solids. Even if starting 
at zero pressure with the same molecular Cmca struc¬ 
ture which CO 2 adopts at moderate pressure, CS 2 need 
not develop at very high pressures the same structures 
as CO 2 . To fill this remaining knowledge gap we con¬ 
ducted an unbiased theoretical crystal structure search 
to explore high pressure structures of CS 2 . This search 
revealed that at high pressures where the carbon coor¬ 
dination is four the /3-cristobalite structural motif is in¬ 
deed superseded by a novel, robust and unsuspected lay¬ 
ered network of tetrahedra. Only while preparing this 
manuscript, fresh high pressure experimental work on 
SiS 2 appeared, 10 reporting the very same structure, des¬ 
ignated there as HP1. Our theoretical search now found 
the tetrahedra-based structure of CS 2 , which is lowest in 
enthalpy in a wide range of pressures from 30 to 170 GPa, 
to be structurally identical with HP1 structure of SiS 2 10 
(see Fig. [I]). Based on this HP1 structure of CS 2 we 
calculated the Raman spectra and pair distribution func¬ 
tions and found it to agree better with experimental data 
than those of /3-cristobalite, the high pressure structure 
of CO 2 We also obtained predictions for the infrared 
(IR) absorption spectra, not yet available experimentally, 
that will hopefully serve in the future to experimentally 
corroborate or discard our predictions. The calculated 
electronic structure of high pressure metallized HPI-CS 2 
GPa moreover indicates a small Fermi level density of 
states, questioning the intrinsic nature of the observed 
superconductivity. 

This paper is organized as follows. In Section [TT] we 
present the details of our structural search and subse¬ 
quent analysis. In Section |III| we present the structural 
results of our search. In Sections IV and[V]we analyze the 
lattice dynamics, Raman spectra and electronic structure 
of the new phase. In the final sectior |VT| we summarize 
the results and draw conclusions. 


II. COMPUTATIONAL METHODS 

The search for low-enthalpy structures of CS 2 was car¬ 
ried out by exploiting an evolutiona ry alg orithm (EA), 
as implemented in the USPEX code. 11 ® The EA was 
run in conjunction with ab-initio electronic structure 



FIG. 1. (Color online) The P2\/c layered crystal structure, 
named HP1 after that reported for SiS 2 p^ viewed along a, 
b and c-axes. Each layer is made up of two pairs of CS4 
tetrahedra. The two tetrahedra in each pair are joined by 
the edge, with S edge atoms shown in red. At « 50 GPa the 
interlayer distance is about 2.1 A. Full structural data are in 
Table [T] 


calculations and relaxations based on standard density 
functional theory (DFT). The Perdew-Burke-Ernzerhof 
(PBEjPSI generalized gradient approximation as imple¬ 
mented in VASP (Vienna Ab-initio Simulation Pack¬ 
age )T^ employ ed th e projector-augmented plane wave 
(PAW) method®® The energy cutoff for the plane-wave 
basis was set to 550 eV to ensure full convergence. The 
Brillouin Zone was sampled by Monkhorst-Pack meshes 
with a resolution of 27 t x 0.05 A -1 . Since at low pressures 
the dispersion forces are important we empl oyed f or en¬ 
thalpy calculations the optB86b-vdW schemd^ISIl based 
on the van der Waals density functional.® Lattice zero- 
point energies were not included and temperature was 
assumed to be zero throughout. The phonon, Raman 
and IR spectra were instead calculated in linear response 
as implemented in the Quantum ESPRESSO suite of 
programs within the LDA approximation. 20 For that, the 
GGA input structures were further relaxed with LDA be¬ 
fore carrying out the phonon calculations. 

Brute-force application of the EA structural search al¬ 
gorithm to a compound such as CS 2 with a positive for¬ 
mation enthalpy should lead by necessity to chemical 
decomposition and outright disappearance of the com¬ 
pound itself. As was said above, even at ambient pres¬ 
sure molecular CS 2 is metastable, meaning it is locally 
stable mechanically, and temporarily surviving due to ex¬ 
ceedingly slow kinetics, while intrinsically unstable on 
thermodynamic grounds. The EA search ignores kinet¬ 
ics and will therefore lead to decomposition, given a large 
enough supercell. In this study, where we wish to find 
and study the high pressure metastable phases of the un¬ 
decomposed compound, a strict CS 2 stoichiometry was 
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TABLE I. Structural data of the new phases. As seen from 
Wyckoff positions, the number of CS 2 formula units/cell is 
Z = 4 for all the phases except HP2 where Z = 12. 
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assumed with an EA cell chosen to contain 12 atoms, i. e. 
four CS 2 molecular units. Even if small, that cell is actu¬ 
ally still large enough to show decomposition and phase 
separation into carbon and sulfur in an unrestricted EA 
search. We therefore artificially prevented decomposition 
by constraining the EA search to avoid the formation 
of C—C and S—S first-neighbor bonds. This allows to 
find the high pressure crystal structures where decom¬ 
position, total or partial, is artificially suppressed. It 
should be noted that while the artificial nature of this 
constraint is such that it may endanger the formation of 
some denser structures at higher pressures, still various 
octahedral structures were also created, but their calcu¬ 
lated enthalpies were considerably higher with respect to 
the tetrahedral ones. Therefore, we will by construction 
limit ourselves to explore metastable structures with co¬ 
ordination larger than four. With this constraint, the 
crystal structure search was performed at 0, 26, 38, 75, 
120, and 170 GPa. The EA initial population (the num¬ 
ber of structures in the starting generation) was set to 
120 , a large number chosen in order to densely sample 
the configuration space of the random search, later re¬ 
duced to 30 in the following generations. 


III. CRYSTAL STRUCTURES, PAIR 
CORRELATIONS AND X-RAY DIFFRACTION 
SPECTRA 

The EA search produced a large variety of struc¬ 
tures. In particular, near P = 0 we reproduced the 
known molecular Cmca structure, validating the scheme. 
Since experimentally the CS 4 tetrahedral structure was 
created at 30 GPa,® we further focused on search for 
tetrahedral structures at pressures above 26 GPa. In 
Fig(2]we present a selection of low-enthalpy phases, that 
is molecular Cmca , and among the non-molecular, a- 
and /3-cristobalite, compared with the tetrahedra-based 
phases HP1 and HP2. To complete the comparison, we 
also added in the a-quartz structure well known from 
Si 02 - While as explained above all these structures are 
metastable against decomposition, those which we will 
describe are at least mechanically stable under struc¬ 
tural optimization and and dynamically stable as shown 
by phonon calculations (see e.g., Fig. [ 6 ]), thus represent¬ 
ing local enthalpy minima. At any given pressure these 
phases compete, so they are all doubly metastable except 
for that with lowest enthalpy. It makes sense to consider 
them all anyway, as that broader picture clarifies the 
relative stability margins - and also because metastable 
phases commonly appear experimentally for kinetic rea¬ 
sons. 

Except for 0 GPa where molecular Cmca was found 
to prevail in agreement with experiment, the most im¬ 
portant structure obtained by the EA search at high 
pressures was the P2\/c (HP1). This monoclinic layered 
structure is quite interesting. We optimized its structure 
at 50 GPa and found the structural parameters as listed 
in Table |T| As shown in Fig. [T] the carbon coordination 
is four, but is quite different from /3-cristobalite. Each 
layer consists of four CS 4 tetrahedra, or more accurately 
two pairs of tetrahedra. Tetrahedra of different pairs 
(different colors in Fig. [l]) share a corner sulfur, but two 
tetrahedra in the same pair (same color) share an edge, 
made up of two sulfurs. Simple as it looks, this structure 
came in at first as a total surprise—we designated it as 
’’shaliabite”, in the lack of an existing name. We found 
subsequently that the very same structure has been very 
recently 10 observed and called HP1 in the phase diagram 
of SiS 2 , at the much lower pressure of 2.8 GPa (actually, 
it had apparently been observed in SiS 2 long time agd^ 
but not resolved). 

In RefflQ] another monoclinic phase with the space 
group P2i/c denoted as HP2 was observed to follow HP1 
at pressure of 3.5 GPa in SiS 2 - Unlike the HP1, which 
is layered and involves four CS 2 units, the HP2 phase 
consists of 12 CS 2 units, arranged in a 3D covalent net¬ 
work. For the sake of completeness we calculated en¬ 
thalpies in CS 2 for optimized HP2, and for the chain-like 
orthorhombic phase with space group Ibam which was de¬ 
noted in Ref{TO] as NP, the structure of SiS 2 at ambient 
pressure. The structural parameters of these theoreti¬ 
cal CS 2 phases (HP1, HP2 and NP) at various relevant 
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FIG. 2. (Color online) Zero temperature enthalpy for some 
optimized structures of CS 2 , plotted relative to /3-cristobalite 
(J42d). No structures with carbon coordinations of three or 
larger than four are considered. Structures HP1 (P2i/c ), 
HP2 ( P2\/c) and NP ( Ibam ) are based on CS4 tetrahedra. 
The relative stability of HP2 between 10-30 GPa is marginal 
and possibly unreliable, a-cristobalite* (P2i2i2i (No. 19). 
) is a high pressure version of regular a-cristobalite (P4i2i2 
(No. 92)) where symmetry is reduced through a rotation of 
the tetrahedral environment of carbon atom. 


pressures are listed in Table. [TJ Among phases not shown 
here and not further discussed is tridymite, also consid¬ 
ered by Ref{7l which we found at 50 GPa to be unstable, 
spontaneously converting into a low-symmetry structure. 

The relative enthalpies calculated for these structures 
are shown in Fig(2] The Cmca molecular phase pre¬ 
vailing at low pressures is quickly supplanted by I42d 
/3-cristobalite at about «10GPa. At the same time, 
however, the i 3 2 1 /c tetrahedra-based structures appear, 
preempting this transformation and replacing the Cmca 
molecular structure already at about 5 GPa. Upon in¬ 
creasing pressure, the two structures HP1 and HP2 re¬ 
main nearly iso-enthalpic up to 30 GPa. Above that pres¬ 
sure the HP1 layered structure clearly prevails, its en- 
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FIG. 3. Calculated X-ray diffraction patterns of the HP1 and 
/3-cristobalite structures of CS 2 at P = 55 GPa, compared with 
experimental data of Ref. [7] 


thalpy remaining lower than that of /3-cristobalite by a 
substantial amount, exceeding the estimated computa¬ 
tional errors of about 10 meV/molecular unit, up un¬ 
til 160 GPa—a regime where, however, carbon coordi¬ 
nations larger than four will take over. 7 ! The same HP1 
structure which we found here appears in SiS 2 at much 
lower pressures than in CS 2 as is natural given the shorter 
bond lengths and smaller compressibility of the carbon 
compound. 

It is interesting to rationalize the finding of edge¬ 
sharing tetrahedra in high pressure CS 2 , which behaves 
similarly to low pressure SiS 2 whose phases consist of 
variously packed tetrahedra. The stability of tetrahedra- 
based phases in a carbon compound such as CO 2 is de¬ 
nied by Pauling’s third rule for ionic crystals which states 
that carbon edge-sharing has a destabilizing effect as it 
brings the positive carbon ions too close together in¬ 
creasing their Coulomb repulsion energy. The question 
is then why this obstacle does not arise in CS 2 . In or¬ 
der to compare the importance of this effect in CO 2 and 
CS 2 we structurally optimized CO 2 in the HP1 structure 
at 50 GPa. In agreement with Pauling, we found for 
this structure a much higher enthalpy of 0.4eV/molecule 
with respect to the stable phase /3-cristobalite, show¬ 
ing that edge-sharing in CO 2 is indeed unfavorable. To 
confirm that this is due to ionicity we calculated Bader 
charge^Smi f or the C, S and O atoms in the HP1 struc¬ 
ture. Strikingly whereas in CO 2 the partial charge of 
carbon is about +2, it turned out to be about —0.5 in 
CS 2 . The bond polarization in CS 2 is not only of small 
magnitude, but opposite to that in CO 2 . This finding ex- 
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FIG. 4. (Color online) Calculated pair correlations G(r) at 
50GPa for /3-cristobalite (green) and layered HP1 (blue) in 
comparison with experimental data at 55 GPa taken from 
Ref. H An empirical Gaussian broadening was applied to 
the G(R) of the respective perfect crystals in order to mimic 
the structural disorder. The Gaussian width was chosen in 
order to get for the C-S bond peak at 1.75 A the same width 
as in experiment. 7 



FIG. 5. (Color online) Calculated phonon spectrum of the 
P2±/c layered structure of CS 2 at 50 GPa. 


plains why edge-sharing in CS 2 does not have the desta¬ 
bilizing effect it has in CO 2 . A second significant differ¬ 
ence between CO 2 and CS 2 originates from the different 
chemistry of oxygen and sulfur which becomes relevant at 
higher pressures where the chalcogen binds two carbons. 
In this configuration, oxygen hybridizes sp , favoring a 
bond angle of about 109°; sulfur instead prefers p orbital 
binding without hybridization, and a bond angle closer 
to 90°. (See a detailed discussion e.g.,in ReflSTl). There¬ 
fore, it is not surprising to find that CO 2 should adopt 
the /3-cristobalite structure where the C-O-C bond an¬ 
gles are 106° and 115° while CS 2 prefers the edge sharing 
tetrahedra with C-S-C bond angle of 90°. 

Fig (3] shows the calculated diffraction pattern of the 
HP1 and /3-cristobalite structures of CS 2 at 50 GPa. 
Comparing the position of the Bragg peaks with those of 
broad peaks of the X-ray structure factor S(Q) in experi¬ 
ment at 55 GPa (Fig.3 in Ref[7j) we find roughly the same 
agreement for both crystalline structures. In order to dis¬ 
cuss and compare more realistically the direct-space pair 
correlation function, we carried out an ab-initio Molecu¬ 
lar Dynamics simulation for HPI-CS 2 and /3-cristobalite 
CS 2 at 300 K and P = 0 GPa, using the VASP code. Figj4] 
shows the results in comparison with the experimentally 
extracted G(3?)f’ 7 Although differences are not dramatic, 
the HP1 pair correlations appear to agree with experi¬ 
ment somewhat better than those for /3-cristobalite. As 
a side result, these simulations also indicated a high level 
of stability and robustness of HPI-CS 2 against thermal 
fluctuations. It is believed that this stability will be im¬ 
portant for later tribocliemical studies which we are plan¬ 
ning. 


IV. PHONONS, RAMAN, AND INFRARED 
ABSORPTION SPECTRA 

In order to ascertain mechanical and dynamical stabil¬ 
ity, and in preparation for spectroscopy, we calculated the 
phonon spectrum of the main HP1 phase that dominates 
the phase diagram of CS 2 for a wide range of pressure. 
As shown in Fig. [5j at P=50 GPa all mode frequencies 
are real and positive, confirming the mechanical stability 
of the structure. Comparison with calculated phonons 
of the /3-cristobalite structure 7 shows that modes of the 
HP1 layered structure are slightly stiffer, although there 
is a fair amount of overall similarity. Phonon calculations 
for HP2, prohibitively expensive because of the large 36- 
atom unit cell, were not attempted, also given the uncer¬ 
tain stability of this phase. 

Based on these phonon calculations, we subsequently 
calculated the Raman and IR absorption spectra of 
the HP1 structure in comparison with /3-cristobalite. 
These spectral calculations require an insulating elec¬ 
tronic structure. Therefore, even if available data are 
mostly at higher pressures, we conducted our spectral cal¬ 
culations at 20 GPa where both HP1 and /3-cristobalite 
still have insulating LDA band structures. Actually, we 
found that LDA at 20 GPa yielded a similar volume to 
that of GGA at 30 GPa, and that the crystal structures 
underlying these LDA spectral calculations are not very 
different from those of our previous 50 GPa structural de¬ 
terminations. As shown in Fig|6] the differences between 
the layered HP1 and /3-cristobalite predicted spectra are 
major. The Raman spectrum of /3-cristobalite has a main 
(twin) peak near 300 cm' 1 , a second main peak near 
400 cm - 1, much weaker features near 600-700 cm - 1, and 
nothing at higher frequencies. The HP1 Raman pre¬ 
dicted spectrum exhibits instead a much larger peak near 
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FIG. 6. (Color online) Predicted Raman spectra of two com¬ 
peting CS 2 structures at 50 GPa, compared with measure¬ 
ments at 50 GPa, 297 K.^ The high frequency secondary peak 
near 800 cm -1 is only present in layered P2\/c and absent 
in /3-cristobalite. Also the low frequency spectrum is better 
reproduced by HP1 than by /3-cristobalite. 



FIG. 7. (Color online) Calculated IR spectra of different CS 2 
structures. Note the stiffer frequencies of layered HP1 com¬ 
pared with /3-cristobalite. Blue arrows in dicate the experi¬ 
mental IR peak positions of molecular CS 2 p 8 * ‘ i ^ 


500 cm -1 and considerable spectral intensity at 700 and 
also 800 cm -1 . Both the 500 and the 800 features agree 
much better with experimental Raman data (50 GPa, 
297 K) PI This proves that HP1 is the dominant phase 
of CS 2 , as opposed to /3-cristobalite, near 50 GPa. 

Besides Raman, high pressure systems should also per¬ 
mit the measurement of IR absorption. We therefore 
calculated IR spectra, which not surprisingly turned out 


to be very different for the layered HP1 and for /3- 
cristobalite. As can be observed in Fig[7J the main pre¬ 
dicted absorption peaks of /3-cristobalite are close to ~ 
200, 350 and 610cm _1 . In the HP1 phase instead, af¬ 
ter a weaker structure between 150 and 200 cm -1 there 
is a large and broad absorption band between 700 and 
800 cm -1 , a range where /3-cristobalite should be IR 
silent. In future data, this unmistakeable difference of 
IR absorption spectra should stand clearly out. Experi¬ 
mental IR data exist apparently only for the low pressure 
molecular structure. The IR peak positions observed for 
the Cmca structure represented by arrows in FigjT] agree 
very well with our calculations. 


V. ELECTRONIC STRUCTURE 


DFT calculations yield the electron band structure of 
all low enthalpy phases of CS 2 . All low pressure struc¬ 
tures are insulating. The DFT-PBE electronic band 
structure of the layered HP1 phase is shown in Fig{8](a) at 
40 GPa, above the critical metallization pressure, where 
the band gap of HPI-CS 2 closes. Metallization of HP1- 
CS 2 occurs at about 30 GPa within the PBE functional, 
an approximation which notoriously underestimates the 
gap and therefore the metallization pressure. We re¬ 
peated the CS 2 calculations using B3LYP hybrid func¬ 
tionals, and found the metallization pressure to increase 
to ~50 GPa, a value now in excellent agreement with 
experiment. It is interesting to note that, owing to its 
initially smaller band gap than CO 2 , CS 2 metallizes read¬ 
ily under pressure after turning from a twofold molecular 
state to a dense fourfold solid, unlike CO 2 which remains 
a wide-gap insulator long after a similar transformation 
into /3-cristobalite. One interesting question arising at 
this point is whether superconductivity is predicted in 
the metallic high pressure state of HPI-CS 2 . Unfortu¬ 
nately our limited resources and the large 12-atom unit 
cell prevented us from calculating the electron-phonon 
interaction A and thus estimating the critical tempera¬ 
ture T c . Nonetheless, a qualitative answer to that ques¬ 
tion, even before any detailed calculations, is suggested 
by direct inspection of the band structure of Fig(8jc) . 
Metallization of HP1 takes place by band overlap, with 
formation of a hole pocket at the k = 0 F point and 
of a corresponding pair of electron pockets, at /c-points 
F and F' near Aq and A but off the A — F — Ao line. 
The electron density of states of Fig{8](a) calculated after 
band overlap is small, not suggestive of a large electron- 
phonon coupling parameter A. As in other cases, our 
observation of band overlap metallization 32 does not of¬ 
fer a strong promise of superconductivity, at least of the 
standard BCS kind. However, the wavevectors T-F and 
T-F', ± (0.3636, 0.0000, 0.4545) are electron-hole nesting 
vectors of HPI-CS 2 near the gap closing pressure around 
50 GPa. It is possible that charge-density-wave or spin- 
density-wave static modulations might appear with that 
periodicity, possibly also accompanied by some related 
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FIG. 8. (Color online) (a) The parallelepiped Brillouin zone 
of the HP1 P2i/c structure! 2 ™ F = (0.3636, 0.0000, 0.4545) 
is a general point in the Brillouin zone, (b) PBE electronic 
structure of HPI-CS 2 at 40GPa, above the insulator-metal 
transition, (c) Bands near the Fermi level showing a single 
hole pocket at T and two electron pockets at F and F' = —F. 


superconducting phase. We are not presently in a posi¬ 
tion to inquire quantitatively into this possibility, which 
would require newer and different approaches. 


VI. CONCLUSIONS 

We presented a theoretical study of high pressure solid 
phases of CS 2 . We discarded the obvious possibility of de¬ 
composition by restricting only to phases devoid of C-C 
and of S-S nearest neighbor bonds, and aiming at uncov¬ 


ering the metastable phases of lowest enthalpy through 
an unbiased genetic structure search. Our main result is 
that, contrary to expectations based on similarity with 
CO 2 where twofold coordinated Cmca eventually turns 
into fourfold coordinated /3-cristobalite, high pressure in 
CS 2 eventually converts Cmca into a different fourfold 
coordinated layered P2 1 /c phase built up of edge-sharing 
pairs of tetrahedra. Recently the very same structure, 
named HP1, has been experimentally reported for SiS 2 at 
much lower pressures? 0 Another tetrahedra-based struc¬ 
ture, HP2, is also stabilized between 8 and 30 GPa, but 
only by a small enthalpy difference with respect to HP1 
comparable with our calculation errors. 

The proposed HP1 structure represents a more plausi¬ 
ble candidate structure for high pressure, fourfold coor¬ 
dinated CS 2 than those considered so far?® Both calcu¬ 
lated pair correlations and Raman spectra agree better 
with existing data than those of /3-cristobalite. It will be 
a challenge for the experimentalists to try to prepare this 
phase in a crystalline state amenable to more accurate 
investigation in the future. To that end we provide sub¬ 
stantial additional predictions, in particular of IR spec¬ 
tra, that should be crucial in identifying the new phase. 
Despite its intrinsic metastability, the HPI-CS 2 structure 
appears exceptionally robust. These qualities make CS 2 
a good candidate system for studies of high pressure sim¬ 
ulated tribochemistry, a project which is presently going 
on in the Trieste group. 

Even though the layered HP1 structure is metallic 
above 50 GPa, it seems unlikely that it should account for 
the experimentally observed superconductivity!- 9 because 
the density of states, and therefore the dimensionless A, 
is likely to remain low after the band-overlap metallic 
state. Although a reasonable hypothesis could be a pos¬ 
sible partial decomposition of CS 2 , with creation of some 
free sulfur filaments or other non stoichiometric products, 
we are not in a position to address that occurrence here. 

We note in closing that at pressures just before metal¬ 
lization, where the band gap closing of insulating HP1- 
CS 2 is indirect, the crystal might develop a narrow 
charge-density-wave or spin-density-wave phase, charac¬ 
terized by a nesting wavevector close to ± FpS Although 
there has been so far no observation of this kind in high 
pressure experiments, this possibility, which we also sug¬ 
gested for MoS2p3! seems worth pursuing. 

Last and perhaps most important, our work provides 
a new link into the high pressure crystal chemistry of 
the archetypal family of IV-VI AB 2 compounds made of 
light elements where previously only some cross similar¬ 
ities between CO 2 , Si02, and SiS 2 were discussed. We 
show that there is some structural kinship of CS 2 not just 
to C0 2 at low pressures, but eventually also to SiS 2 at 
higher pressures. 
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